A general treatment of geometric phases and dynamical invariants 
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Based only on the parallel transport condition, we present a general method to compute Abelian 
or non-Abelian geometric phases acquired by the basis states of pure or mixed density operators, 
which also holds for nonadiabatic and noncyclic evolution. Two interesting features of the non- 
Abelian geometric phase obtained by our method stand out: i) it is a generalization of Wilczek and 
Zee's non-Abelian holonomy, in that it describes nonadiabatic evolution where the basis states are 
parallelly transported between distinct degenerate subspaces, and ii) the non-Abelian character of 
our geometric phase relies on the transitional evolution of the basis states, even in the nondegenerate 
case. We apply our formalism to a two-level system evolving nonadiabatically under spontaneous 
decay to emphasize the non-Abelian nature of the geometric phase induced by the reservoir. We 
also show, through the generalized invariant theory, that our general approach encompasses previous 
results in the literature. 
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Introduction. The concept of geometric phase (GP) 
was transposed to the domain of quantum systems un- 
dergoing cycUc adiabatic evolution by M. Berry after 
having been introduced by Pancharatnam in connection 
with interference of light waves with distinct polariza- 
tions 1^. After Berry's discovery, Aharonov and Anan- 
dan [3| removed the need for adiabatic evolution and 
Samuel and Bhandari Q extended the concept of GP 
to noncyclic and nonunitary evolutions, introducing the 
notion of geodesic closure in the projective Hilbert space. 
In recent years, the possibility of achieving fault-tolerant 
quantum computation [1, Q invoked the controver- 
sial subject of GPs for open quantum systems, where 
the dynamic is generally nonunitary. In this case the 
GPs have been defined by different methods: using phe- 
nomenological approaches [8|, stochastic fields [9|, the 
usual master equation and quantum jumps tech- 
niques, apart from state purification mean values of 
distributions Il3|, interferometric 14 1, and superoperator 
methods [il,[li,[l3. 

Parallel to these advances in the understanding of GPs, 
the dynamical invariants (DIs) proposed by Lewis and 
Riesenfeld to handle time-dependent Hamiltonians, 
have been applied to a number of problems, inclu ding 
recent advances in cavity quantum electrodynamics [l9| 
and Bose-Einstein condensates [23|. Based on the DIs, 
Morales 21 1 and Mizrahi introduced, independently, 
a convenient way to compute the evolution of the GPs, 
which was used to suggest an interferometric experiment 
to measure GPs induced by a Stark shift in cavity quan- 
tum electrodynamics [23 1. 

In this letter, we rely only on the parallel transport 
condition to obtain a general formal expression to com- 
pute the GPs acquired by the basis states of a density 
operator under unitary or nonunitary, adiabatic or nona- 



diabatic, and cyclic or noncyclic evolutions. The DIs ap- 
plied to the density matrix [^j are required to account 
for the time-evolution of the GPs in a general scenario 
of open quantum systems. Apart from showing that our 
method reproduces previous results in literature, we use 
it to compute the GP acquired by a two-level system 
under spontaneous decay. We show unambiguously that 
this GP turns out to be non-Abelian even in nondegener- 
ate system. Thus, it is clear from this general approach 
that the non-Abelian character of the GP is wholly asso- 
ciated with the transitional dynamics of the basis states. 
Moreover, our treatment generalizes Wilczek and Zee's 
[25! non-Abelian holonomy, in that it describes nonadi- 
abatic evolutions in which the basis states are parallelly 
transported between distinct degenerate subspaces. 

Parallel transport condition. To compute the GPs 
acquired by the time-dependent states {|A,a;t)} of 
an orthonormal basis, where A is g-fold degenerate 
(a = l,2,...,g), we first introduce the condition 
for parallel transport of a state vector, namely that 
II (A', a'; t| dt |A, a; — [2^. The evolution operator 
responsible for the parallel transport during the time 
interval t, \\,a;t)^^ = Vl^i(i) |A, a; 0), may be written as 

= EA,A'Ea,a'W,)rA'' W|A,a;i)(A',a';0|. Sub- 
stituting V„{t) in the parallel transport condition, we 
find the equation for its coefficients in the matrix form 

V„{t) = iA{t)V„{t), whose solution is given by 



V„{t) = Texp 



drAir] 



(1) 



where ^,(0) = 1, the elements of the non-Abelian con- 
nection A{t) are yl^'^, (t) — « (A, a; r| c?^ |A', a'; r), and 
T is the time-ordering operator. The unitary matrix 
V„(t) {A(t) being Hermitian) accounts for a cyclic and 
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non-Abelian GP. As we are also concerned with non- 
cyclic evolutions, we must account for the relative phase 
acquired by the parallelly transported state \X,a;t)^^ 
with respect to its starting point |A,a;0), given by arg 
{(A,a;0| A,a;t)„} = arg {{X,a;0\W{t,0)V„{t) \X,a;0)}, 
with the elements of the overlap matrix W reading 
W^;y{t,0) = (A,a;0| X',a';t). We note that, to trans- 
port a subspace without rotating it locally, it must re- 
main parallel to itself during the time interval from t 
to t + St, i.e., W^',{t + St,t) = {X,a;t\ X\a';t + St) ~ 
S\,\'Sa,a'- Consequently, a necessary and sufficient con- 
dition stipulates that ^ (<) = 0, and hence A{t) — 0. 

As observed in Ref. [27|, in the general case, where the 
overlap matrix W{t, 0) is restricted to an incomplete sub- 
space ^X^A a l'^' "^i ^) ("^j °' ^1 7^ ' it must be decomposed 
in the polar form W = RU, R being a positive-definite 
(Hermitian) matrix (det R > 0) and U a unitary matrix 



(detC/ 



such that e"'' = det W/ det R. In the case 



where i? is a positive-semidefinite matrix (det R — 0), 
the Moore-Penrose pseudoinverse must be used to evalu- 
ate the matrix U. Then, the non-Abelian, nonadiabatic 
and noncyclic GP acquired by the basis states {\X,a;t)} 
turns out to be 



O{t,0)^U{t,0)V,{t), 



(2) 



where U is the unitary part of W. Although Eq. ([2]) 
has already been published by Mostafazadeh [2^ and 
Kult et al. [l^l, in both cases the analyzed non-Abelian 
geometric phases arise from the usual degeneracy, i.e., 
from transitions inside degenerate subspaces. In a more 
general scenario, we have extended the analyses in Refs. 
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28| to cases of transitional dynamics which connect 
different degenerate and/or nondcgcncrate subspaces. As 
discussed below, we found that transitional dynamics, 
even between nondegenerate states, is the source of a 
non-Abelian geometric phase. Moreover, when consider- 
ing the case of transitional dynamics connecting different 
degenerate subspaces, our analysis reveals a generaliza- 
tion of Wilczek and Zee's non-Abelian holonomy, since it 
describes nonadiabatic evolutions where the basis states 
are parallelly transported from one degenerate subspace 
to another. The definition of parallel transport used here 
appear similar to that used in Ref. [29j , however, our ge- 
ometric quantity ^ is different from that one defined in 
[29! . It is important to note that the treatment presented 
here is general and includes several kinds of system dy- 
namics. 

The general expression for O(i,0) can be simplified 
in three particular cases: transitional (t) dynamics in 
a degenerate (d) subspace (O*''i(<,0) = 0)K,(i)]*''^) 
and transitional 0) = [[/(t, 0)V:,(t)]*''''^) and non- 

transitional (nt) (O"''"'i(t,0) [C/(t,0)V;,(t)]"*'"'^) dy- 
namics in a nondegenerate (nd) subspace. The three 
different cases mentioned above, covering all possible evo- 



lutions of the basis states, will be explored below in the 
context of the DIs for a general Lindblad evolution. 

Gauge-invariance of the GP. Under the gauge trans- 
formation |A, a;t)' = X^i/ b -^^I'J a(*) 1^' where M 
is a unitary matrix, the operator V„{t) becomes 
= MHt)V,{t)M{0), with A'{t) = M'i{t)A{t)M(t) + 

iM^{t)M{t). In its turn, C/(t, 0) becomes f/'(t,0) = 
M\Q)U{t,0)M{t). From these transformed expressions, 
it follows that C(i,0) = Mt(0)O(i, O)M(O), making the 
trace and the eigenvalues of C(t, 0) observable gauge in- 
variants. 

Summarizing our method, we first choose the time- 
dependent basis to expand the state of the system p{t). 
Next, we analyze the evolution of the basis states to ver- 
ify if there are transitions or not between them. With 
these informations we build, through Eq. (2), the op- 
erators O'''i(t,0) = [C/(t,0)K,(i)]*''', or C'"d(t,0) = 
[C/(t,0)K,(i)]''"^ or O"''"d(i,0) = [[/(t,0)K,W]"''"'^, 
given the geometric phases acquired by the basis states 
under a general evolution. 

GP and the DIs. To illustrate the appearance of the 
GPs for the different types of evolution discussed above, 
it is convenient to assume the basis states {|A, a;i}} to 
be the eigenstates of an Hermitian DI, I{t) \ X,a;t) = 
A(i)|A, a;i}, since the invariant method is applicable 
for adiabatic and nonadiabatic processes under uni- 
tary and nonunitary evolutions. Moreover, the infor- 
mation about open quantum dynamics is naturally car- 
ried through the eigenstates of the invariant operators. 
In fact, the DIs 18| associate with any density oper- 
ator p{t) the time-conserved quantity I{t), called in- 
variant, which satisfies the condition d{I{t)) / dt = 
[il ]. From the Lindblad form for the density op- 
erator of an open quantum system described by the 
Hamiltonian Ho{t), given by dtp{t) — —i[HQ{t),p{t)] + 

%j{t) { [r„ p(t)r]] + [r,p(t), r]] }, the evolution 

of the invariant operator following from {I{t)) = 
Tr [I{t)p{t)] becomes 

dj{t) = [Ho{t),I{t)] + 7,,(t) {r] [r„/(0] 

-f r]]r.}, (3) 

where F^ are the Lindblad operators coming from the 
action of an environment and 7ij(t) are the coupling 
strengths. We observe that the Lindblad equation pre- 
sented above applies for the case where the time scale re- 
lated to the rate of change of the time-dependent Hamil- 
tonian Ha (t) is much slower than the time scales charac- 
terizing the reservoir (30| . 

Formal solution of the master equation and its connec- 
tion with GPs. Expanding p{t) in the eigenstate basis of 
the invariant, we obtain through the master equation in 
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the Lindblad form, the coupled differential equations 

A, a 

A, a 

+25^7., E ^Ai'iK'i.^*^ (4) 

A, A', a, a' 

where we have defined the time-dependent ma- 
trix elements H^"^^ = — {ii,b;t\Ho\X,a;t), A''^^ ~ 
{li,b;t\idt\X,a;t), d|^''^ = X]ij7y (m, ^; *l r]r, |A, a; i), 

and AJ^'*^.^ = {fi,b;t\Ti\X,a;t). It can be verified that 
the formal integration of Eq. ([4]) leads to 

c{t)=Tcxp{i ( dT[H + A + iD]{-) + 2y2. ■ [ dr 
[ Jo *J Jo 

X 7yA,(-)A] dT[H + A-iD]jc{0), (5) 

with (•) indicating the side where the matrices H+AzLiD, 
Ai and A], are supposed to act on c(0), in the time- 
ordered expansion of the r.h.s. of ([5]). Differently from 
the original procedure by Berry [1] (within the context of 
an adiabatic and cyclic evolution of a pure state) , it is no 
longer clear, as first noted in Ref. [2^, how to extract the 
GP from the dynamics of the probability amplitudes (Eq. 
([5|)), describing the nonadiabatic evolution of an open 
system, since in general A does not commute with H, D, 
and Ai. Some exceptions correspond to the s pec ific cases 
of open quantum systems discussed in Refs. [l5|, [l7| and 
filtering evolutions [sij. That is why we have presented a 
formal definition of the GP in Eq. ([2]), instead of trying 
to obtain it, as usual 0, 0, S, [32] , from the dynamics of 
the probability amplitudes. For some particular situations 
of dissipative-free dynamics, a connection can be estab- 
lished between the GPs emerging from our definition and 
the dynamics of the probability amplitudes coefficients, 
as shown below. 

Dissipative-free dynamics (degenerate case). For 
dissipative-free dynamics [7ij(t) = in Eq. ^] the 
solution ([5]), inside the degenerate subspaces fi and 
reduces to 

xTe-^^o'^(^"' (6) 



where c,. 



and A 



posed of the elements 



b.a 



Ref. [33], the assumption that [Hi 



represent matrices com- 

Tb,a Ab,a 

= 

even for a nonadiabatic evolution, resulted in the 

rt 

expressions 



and AJ;-,;. 

At, At) ^Ai.Ai] 



Texp 



Texp 



ever, as in general 



Jp dTi?^_^(r) for dynamic and 

for cyclic geometric phases. How- 
Hfi^fi, Afj_,fj] 7^ for transitional 



and nonadiabatic dynamics in a degenerate subspace 
/i, we must return to the formal expression ^ with 

K'-'^W - E^Ea,a' Wi);:' (*) to Obtain 

the GPs for cyclic {U — 1) and noncyclic {U ^ 1) evolu- 
tions. For the cyclic case, Eq. ([2|) with T/*''^ 



{t) leads ex- 
actly to the expression T exp i drAf^^i^^r) [3l| for the 

GP, whereas for the noncyclic case the matrix U must be 
taken into account (instead of matrix W as in Ref. [iSl), 
with the elements of the non-Abelian connection being 
given by A°;^'(r) = i {fi, a; r| dr \fi, a'; t) . 

For transitional and adiabatic dynamics in a degen- 
erate subspace, the adiabatic evolution V^''^{t) is built 
up from the condition dtl — [l^l, which implies that 
Hq and I assume the same basis {|^, a; t)}. In this case, 
the non-Abelian dynamic phase reduces to the Abelian 
one, since H^'^ = E^^Sab, where E^j^ is the eigenen- 
ergy associated with the state |/z,a;i) and, consequently, 
[Hfi^fi, Af^,fj] ~ 0. The cyclic GP thus remains the non- 
Abelian holonomy Texp i dTA{T) computed in the 

Hamiltonian eigenstates, obtained by Wilczek and Zee 
[25! . For a noncyclic evolution, the nonidentity matrix U 
is responsible for yielding the result in Ref. [27| . 

Dissipative-free dynamics (nondegenerate case). 
For the nondegenerate case, we easily verify that 
the solution of ([3]), corresponding to a partic- 
ular case of Eq. is given by c^,^'(i) — 



c,i,,i' (0) exp 

I, _ I,' — H I, _ I, Hi,' 11/ . H I, _ I, — 



AA 



so that 
the Abelian 



/^,/^ ^^M' ,At' ' 

the dynamic and 



and A 



cyclic 



-)]}, with 

{fi;t\Ho \fr,t), 
= {l.r,t\idt\fJ.;t), 
GPs reduce to 



expressions exp 



exp 



and 
and 



Note that the dynamic 
geometric phases associated with the element p^,^' is 
simply the difference between the phases acquired by 
the states \fJ.',t) and even for the nonadiabatic 

case, as can be verified from the DIs (3]. This strik- 
ing feature of the DIs enable nonadiabatic evolutions 
with non-transitional eigenstates in the nondegenerate 
case. In this connection, the operators [/"''"^(t, 0) and 
^„t,nd(^) = Ea (K)aa (t) l^;<> (A;0| in Eq. m describe 
the GP for noncyclic evolutions whose expression is 
(A;0 |A;i)exp(^-/(,*dr(A;r|d^|A;T)) For cyclic 

evolutions this GP reduces to that obtained in Refs. 
21, 22|. Evidently, for the adiabatic case (dtl — 0) the 
above expression for the noncyclic GP gives the result 
in Ref. [33|, while in the cyclic case it reproduces the 
original Berry phase 

Application - An unstable two-level system. A partic- 
ular form of the operator V„{t), given by [V^i(t)]''"'^ = 
Ea A' (^i)aa' (*) 1^! (-^'i 0|, occurs when we have a 
nonadiabatic evolution with a transition between non- 
degenerate states of the basis {|A;i)}. This is the case 
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in a quantum system undergoing a unitary evolution, as 
in Ref. (26| . or a nonunitary evolution where, as shown 
below, the noise injection gives rise to a non-Abelian GP. 

Employing the above method to analyze the role of 
dissipation in the evolution of the GP, we consider a non- 
degenerate two-level system, with transition frequency 
luq between the ground (g) and excited (e) states, under 
spontaneous decay at OK. The dynamics of this system 
is described by the solution to the master equation for 
p{t) with Hq = wo<Tz/2, 7ii = 7/2, and Fi = (t_. For 
the invariant operator, we assume 



lit) 



E„ Xap{t)(Jal3, 
o.J5—g^e 



(7) 



where aap — \a) {I3\ and the coefficients Xapit), which 
are solutions of Eq. ([3]) , satisfy the relations 



Xgg{t) = -ro cos 6*0, 
Xeeit) = (2e^*-l)rocos0o, 
Xeg (t) = ro sin do e 



Xgeit) = rosin6'oe 



tt/2-i(uiot+M 
7t/2+i(cJot+0o) 



(8a) 
(8b) 
(8c) 
(8d) 



with 9o (00 ) being the polar (azimuthal) angle of 
the initial pure (ro — 1) or mixed (ro < 1) state 
p(0) in the Bloch sphere. The eigenvalues and eigen- 
vectors of the invariant (O are given by X±{t) = 
-ro \{1 - eT*) cos 0o T e^*/^ ^1 - (1 -e'r*)cos2 0o] , 

\±-t) = ±7V(t){/(O|^)±rosin0ocT'(-°*+^")|;)}, 

with f{t) = -ro [eT*/^ cos go - ^1 - (1 - e''*) cos^ 60 



-ro 

and N^t)^ {p{t)+r^sm^0oy\ 

In order to find the coefficients (Ki)aa' (*) '^^ 

n t,nd 



t,nd 



we first solve the equation 



t,nd 



where the elements of the non-Abelian connection are 
given by Aki = i {k;t\dt\l;t), with k,l = +,-. 
To this end, we move to the rotating frame R — 
exp [77 (e~*^ (T_ — e"» (t+) /2] , as in Ref. [2(|, obtaining 
the coupled differential equations 



2A^^ro sint 



oJo/sinO + / cos 8 



2N'^ro sin 9o {1-1^0^0 sin 9o 



cot 77 



-cJo/cosG + / sine | , 



(9a) 



(9b) 



where O = LUot + (j)o — C,- Under the initial condition 
■"'/T' ^ ^0 ^ (which is necessary to obtain an approx- 
imated analytical solution), the assumption of a typical 
weak system- reservoir coupling, 7/^0 ^ li a-nd assum- 
ing a time evolution around t 2± 27r/a;o, we obtain r) ~ 0, 

such that ri{t) ~ arccot |cot6'o 1 + ^^I'-cosflo^'' } ^^'^ 
C{t) ^ ujot + (j)o - -f cos 6'o/2cjo- 

^From the above result, we obtain the non- 
cyclic, nonadiabatic, and non-Abelian GP 



associated with the decaying two-level sys- 
tem O*'"'i(i,0) = L/(i,0)i?t(t)e*"(*)'^- i?(0), 
where n{t) ~ woi [1 -I- 7S'(6lo)t/2], 5(6*0) 
-cos6lo (1/2 - cos6'o + 3cos2 6*0/8) /(I -cos 6*0), and 



C/(i,0) 



UoD U'h 



(10) 



with the on- and off-diagonal elements of 
the overlap matrix given, respectively, by 
Ud ^ N sin (6lo/2) [/ + 2ro e-^'^o* cos^ (6*0/2)] and 
UoD = N e-">"> cos (6I0/2) [/ - 2ro e"'""* sin^ (6lo/2)] . At 
this point, we stress that the decay process introduced 
by the reservoir leads to a transitional dynamics of 
the DI basis states, bringing about a not fault-tolerant 
non-Abelian GP. In the particular case 7 = (nonde- 
generate dissipative-free dynamics), the non-transitional 
evolution during the time interval 2t:/luo leads to the 
cyclic GPs C'^^*'"'^(27r/a;o, 0) = ±7r(l — cos^o) associated 
with the eigenstates |±; t), as obtained in Refs. ([H, H^]). 

Whereas some works are concerned with a formal def- 
inition of the GPs for open quantum systems [13, [3, 
14, 3, 3 12 1, most of them restrict themselves to com- 
pute corrections to this phase coming from the reservoir 
E S E EH- From the formal approach presented in 
this work we additionally verify that the reservoir may 
even change the nature of the holonomy, from an Abelian 
(7 = 0) to a non-Abelian (7 7^ 0) one, as it emerges from 
the above application. 

Summarizing, we have presented a general formalism 
to compute GPs, starting only from the parallel trans- 
port condition. These GPs transform covariantly and 
the approach is applicable to a general scenario, including 
adiabatic or nonadiabatic, cyclic or noncyclic, and tran- 
sitional or non-transitional evolutions of pure or mixed 
states. Although we have used the DIs to compute the 
GPs acquired by the basis states of the invariant, the for- 
malism is applicable to any time-dependent basis states. 
Besides reproducing well-known results established in 
the literature, our formal definition reveals two strik- 
ing features of the GP: it g eneralizes Wilczek and Zee's 



non-Abelian holonomy [25|, in describing nonadiabatic 



evolutions where the basis states are parallelly trans- 
ported between distinct degenerate subspaces; secondly, 
our method demonstrates clearly that the non-Abelian 
character of the GP arises from transitional dynamics, 
even in nondegenerate case. We have shown that the 
nonadiabatic evolution of an open two-level quantum sys- 
tem introduces a non-Abelian holonomy. Both of these 
features not only deepen our understanding of GPs, but 
also offer the possibility of investigating how to use the 
non-Abelian holonomy acquired by transitional dynam- 
ics between nondegenerate states to perform geometric 
quantum computation. 
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